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ABSTRACT 


Computer programs are designed to calculate and display beam patterns in both 
two and three dimensions. Graphical presentation and evaluation in three-dimensions 
are difficult and important problems. 

Five computer models are presented and used in investigating the directional and 
frequency properties of shaded and phased doublet, triplet and quadruplet arrays. 
Comparisons associated with parameters such as wavelength, inter-source distance, 


source strength, and phase difference are examined. 
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I. INTRODUCTION 


Not only is the physical interpretation of the radiation of the acoustic waves in 

more than one dimension complex, but the effects of varying frequency and shading 
are hard to visualize when considering directionality and frequency properties. As 
Wood implied in his book, A Textbook of Sound (Ref. 1: p.449], the directional factor 
is very important in many acoustic applications such as speech, hearing and array 
design. 
The simplest array one can consider is the beam pattern of a dipole radiating in two 
dimensions. References to the case of the simple dipole can be found in virtually every 
textbook in acoustics. For a fairly recent example, see Malecki’s book [Ref. 2: p.199]. 
Currently, the mathematical modeling of the simple acoustic arrays has been enhanced 
with the advent of the sophisticated computer software. Specifically, Ziomek 
(Ret. 3: p.138] and, Embleton and Thiessen (Ref. 4: p.1124] have investigated N point 
source arrays, considering both amplitude and phase weighting. Although the 
mathematics and computer hardware exist, the visualization of the three-dimensional 
aspects of these arravs 1s difficult. By use of the computer graphics software the ability 
_to interpret the significance of shading and phasing can be fully developed (Refs. 5,6]. 

To date, visualization of the three-dimensional beam pattern has been needed to 

understand the directional factor. Numerical computer output was essential to draw 
three-dimensional graphs. By changing two parameters (in three-dimensions, those 
could be x and v or 8 and Q, etc.), different generated values of the beam pattern, F, 
are obtained. 
One computer program to visualize the beam pattern in three-dimensions (P3D) and 
four computer programs to read values in detail (2DICIR, 2DPCIR, 2DIFIX, 
2DPFIX) were written. These computer models predict the directional factor and 
overall beam pattern with or without differences in the source strengths, inter-source 
distances, and phases of the elements of the arrays. 

Typical geometry for the beam pattern prediction is shown in Figure 2.3 on page 
16, 


The following definitions apply and will be used throughout: 


I = source strength of the first source (source 1) 


A = source strength of source 2 


B = source strength of source 3 
C = source strength of source 4 
D, = distance from the geometrical center to source | 
D, = distance from the geometrical center to source 2 
D, = distance from the geometrical center to source 3 
D, = distance from the geometrical center to source 4 


@, = phase difference between source 2 and source | 

; = phase difference between source 3 and source I 

@, = phase difference between source 4 and source | 

§ = angle measured from the z axis toward the x-y plane 

( = angle measured counterclockwise in the x-y plane starting from the x axis 


k = wave number 


w = angular frequency 


R = distance from the geometrical center to a certain point in space 
H = directional factor 


F = beam pattern = 20log,,|H| or 10log, jH? expressed in dB 


All distances are measured from the origin of the coordinate system to the sources 


along the x and y axes. Phase differences are based upon the phase of the first source, 


which is assumed to be zero. 


The goal of this thesis is to develop a convenient methodology for the 


investigation of the functional. dependency of the directional factor on the source 


amplitudes (A, B and C), distances from the origin (D,, D5, D, and D,), phase 


differences (5, @, and @,), and 9 and ® as demonstrated by the five programs: P3D, 


ZDICIR, 2DPCIR, 2DIRING 2DPRIS 
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Il. THEORY 


A spatial configuration of simple sources, either discrete or distributed, each with 
its own complex source strength, can be used to represent the quadruplet source case. 
The pressure at a field point is the sum of pressures produced by the individual sources 
[Ref 7: p.169]. Before we develop the quadruplet source case, the acoustic doublet 
case needs to be described to provide a background for the more complicated 


quadruplet case. 


A. DOUBLET SOURCES 
A doublet source consists of two point sources of strengths of 1 and A, separated 
from the origin, by distances D, and D> and vibrating at the same frequency but with 


a phase difference o. 





Figure 2.1 Geometry of doublet. 


From the geometry shown in Figure 2.1, we have 


cos® = 6,/D,, (2.1) 
sind = ¢,/D,, (2.2) 


i = V(R - ene 5 6, 


= A(R : D, sin@)? + (D,cos@)’, (2.3) 
— = J(1- D,sin0/R)* + (6,cos0/R)’. (2.4) 
sind = ¢,/D,, (2.5) 
cos0 = 05, D (2.6) 


m= J(R + g,)° + 6,° 





= /(R + D,sin@)* + (D,cos®)’, Qa. 
a J(1 + D,sin®/R)? + (6,cos®/R)?. (2.8) 


re 


The pressure at point (R,®) due to source | is 


l ; 
r 
] 


where w is the angular frequency. That due to source 2 is 


A : 
Paci ee (2.10) 
2 


The total acoustic pressure at point (R,@) is then the sum of P, and P., 


p = —— ei(wt-kR), R jk(R-1,) + AR Sjtk(r,-R) + 1) (2.11) 
R t r, 


where R is the distance from the field point to the origin. From (2.4) and (2.8), we get 


k(R - t,) = kR(1 -1,/R) 


Ke 


= kR{I - ¥(1 - D,sin@/R)? + (D,cos0/R)* J, (2.12) 
k(t, - R) = kR(r,/R - 1) 


= kR[V(D,sin®/R + 1)? + (D,cos@/R)* - 1]. (2.13) 


Therefore, equation (2.11) can be rewritten as 


pe — oj(wt - kR) AR jkRUI - Y(1 - D,sin@/R)? + (D,cos0/R)? } 
R r 
] 


, OR ikRE/(D,sin®/R + 1)? + (D,cos0/R)?- 1] + 9] \ (2.14) 


7 
In the most frequently encountered cases, the observation of the pressure is made 
at distances greater than the separation of the sources. Therefore, it will be useful to 
derive the far field (R> >(D, +D,)) pressure equation from the above equation (2.14). 
From the geometry shown in Figure 2.2, we have 
sin@ = €,/D, = £,/D,. (2715) 


a 


Therefore, (R - r,) and (r, - R) can be replaced with D,sin® and D,sin@. Also we have 
—_ ~ —_ ~ 1 (2.16) 


so that the pressure equation becomes 


p = —el(Wt- KR) ikD,sin® 4 4 ¢i(kD,sin6 + 9) (2.17) 
| R 


IP| = P(R)H(@) (2.18) 


where P(R) = (1+A)/R. 
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H(8) = See e(kD,sin®) + aei(kD,sin8 + 9) | (2.19) 


This directivity factor, H(®) can be written as 


_ 2 
H(8) = Tita) ee + Acos(kD,sin8 + 9) 


+ [sin(kD, sin®) - Asin(kD,sin8 + @)]?} 7. (2.20) 





Figure 2.2 Geometry of doublet (far field). 
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B. QUADRUPLET SOURCES 
Similarly, a quadruplet source case can be developed with four sources of 


strengths 1, A, B and C, separated from the origin, by distances eh D,, D, and D 4 
and vibrating at the same frequency with phase differences ®,, @, and @, (Figure 2.3). 





Figure 2.3 Geometry of quadruplet. 


From the geometry shown in Figure 2.3, we have 


= RsinOcos®, Qo 
= Rsin@sin®, (2.22) 
= Reos8, (2.23) 
R= /x*+y*+2? (Oy 


where x, y, Z are rectangular coordinates and R is the distance from the field point to 
the origin. 
And then 


fila JA(y + Ds +x? + 2°, (225) 


By Taylor’s approximation 


Ap oe Dy ere Re (2.26) 
ee Dy ee D, sinOsin®, (22m) 
f> = AyD.) x? eee (2.28) 
R-, = D,y/R = D,sinOsin® (2.29) 


where R is much greater than D,, and D,. 
Using the same approach 


ty = J(x-D,) = yy Seay See ee 


= RV1 - (2D,x/R* - D,*/R*). (2.30) 
r, ~ Ril - (1/2)(2D,x/R? - D,?/R*)] = R-D3x/R + D,7/2R, (2.31) 
R-r, = D,x/R = D,sinOcos® (2.3725 
tr, = V(x + Diy yy + oe Ree, (2.33) 
r,-R = DysinOcos® (2.34) 


where R is much greater than D,, and D,. 
Therefore, the pressure at (R, 9, ) is 


p = + ,i(wt-kR-kD,sinOsin®) , A oj(wt-kR + kD,sinOsin® - @,) 
R R 
sie B el(wt-kR + kD,sin@cos® - @,) 


ais C e]( wt -kR- kD, sinOcos® = M4) (235) 


The directivity factor, 


H(0,0) = ! | ec(kD,sin@sin®) 4, 4 ,i(kD,sin@sin® - ¢,) 
(1L+A+B+C) 
+ Bel(kD3sinBcos® - 3) + Ceri(kDysinBcosD + ,) | (2.36) 


Expanding the complex terms yields 


l 
H(8,@) = ta ae + Acos(kD,sinOsin® - 9.) 


+ Bcos(kD,sinOcos® - @,) + Ccos(kD,sinOcos® + ,)I° 
+ [sin(kD, sinOsin®) - Asin(kD,sin@sin® - @,) 


- Bsin(kD,sin@cos® - @,) + Csin(kD,sinOcos® + ,)1°} 12 (ese 


This 1s the expression used for the computer program. 


C. GRAPHICAL PRESENTATION OF THE DIRECTIVITY 
1. Two-Dimensional Rectangular and Polar Plots 
From equation (2.37), H (0@,@) versus 0, and H (0,@) versus @ can be plotted 
if either ® or @ respectively- is given constant value. In a rectangular plot, the 
directivity, F, which is equal to 20log,,|H (6,@) | is plotted on one axis and 8 or M on 
the other. An example of this is seen in Figure 3.3. In a polar plot, the directional 
factor or beam pattern is a magnitude of the radial distance from the center of the plot 
while the appropriate angle is measured clockwise around the origin. An example of 
this 1s seen in Figure 3.4. 
2. Three-Dimensional Cartesian Plot 
From equation (2.37), H can be calculated and plotted as a function of both 8 
and ®. This provides the three-dimensional surface with its elevation equal to 
20log, glH (8,@)| for each 8-@ pair. An example of this is Figure 3.7. This form of 
presentation is difficult to interpret but can be a very useful analytical tool. Further 
explanations are given in the next chapter. 
3. Contour Plot 
A contour plot gives lines of constant H (0,®). As an example, let kD, = 
See KD = KD, — KD Oo, — 0,0, = 0, = mand A = B = C = lI. Equation 
(2.37) then can be simplified to 
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l 
H(8,M) = =a lcos(kDsin@sin®) - cos(kDsin@cos®)|. (2.38) 


By setting H(6,M) = Q, the desired value, and sin@sin® = c,, and sinBcosD = C5 


where [c,| = ICI, we get 
Icos(kDc,) - cos(kDe,)| = 2Q. (2.39) 


This transcendental equation is then solved for c, and c,. Now, with c, and c, 


determined, we need to solve the second transcendental equation 


o 
tan@® = 





(2.40) 
Cate 


which was derived by dividing c, by c,, for ® and then, for 0 using the 0-© 
relationship. However, many of c,-c, pairs are possible solutions to the first 
transcendental equation. Therefore, equation (2.40) must be solved for all the different 
C,)-C, pairs by repeating the above procedure. This would result in a contour plot such 
as Figure 3.10. This analytical approach was not pursued further in this research 
because other algorithms for extracting constant values of F for varying 8 and ® had 
already been obtained and successfully applied. However, the utility of the above 
transcendental equation approach may well bear further investigation, and we will 
pursue a little case-by-case analysis for some cases for which the solution of equations 


(2.39) and (2.40) is relatively straightforward. 
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Ill. PROGRAM DEVELOPMENT 


In analyzing the directional and frequency properties of shaded and phased 
simple arrays several types of display methods were needed. At the first stage of 
developing the computer program, a HEWLETT-PACKARD 86B computer and a 
HP7090A MEASUREMENT PLOTTING SYSTEM were used. By using this 
equipment we were able to get two-dimensional plots as shown in Figures 3.1 and 3.2. 
However, it took much time in executing the program and also we could not do three- 


dimensional plottings. 
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Figure 3.1 2-D rectangular plot (HP 86B). 
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Figure 3.2 2-D circular plot (Hiv seu) 


To get a complete understanding of the beam pattern, three-dimensional displays 
viewed from various angles were essential. Consequently, alternative methods of 
display were programmed in FORTRAN and executed on the IBM 3033. There are 
several reasons why this proved to be the most desirable method. Primarily, a three- 
dimensional graphics software package (DISSPLA) is available for use on the IBM 
3033. Secondly, FORTRAN generated the numerical output with greatly increased 
efficiency when using the EXTENDED FORTRAN COMPILER [Ref. 8: p.18]. 
Finally, the IBM 3033 has many more output devices such as VERSTEC, TEK618 and 
SHERPA to name a few. Therefore, the programs used in the HEWLETT- 
PACKARD 86B were rewritten in FORTRAN. The translation from BASIC to 
FORTRAN required some changes in variable names along with the use of some 
convenient function calls. For reference, some of the original BASIC programs are 


attached in Appendix A. 


THETA Vs. DIRECTIVITY 


a DEGREE 
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eS 
= 
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Figure 3.3 2-D rectangular plot (DISSPLA). 


In the trial of analyzing only with numerical outputs, BASIC had become 
ineffective in recognizing details due to the very large number of computed values. At 
this stage the graphic package DISSPLA, available on the IBM 3033 computer, was 
used in conjunction with these programs. The DISSPLA package is a library of 
FORTRAN subroutines that facilitate data plotting. It is also device independent, 
which means that it does not rely upon features particular to any type of graphical 
device and has its own extensive symbol and character generation routines. By being 
device independent, programs of the HEWLETT-PACKARD 86B could be translated 
without difficulty. Therefore, DISSPLA was used to produce both two-dimensional 
and three-dimensional graphs. The data used to compute graphs could be generated 
within the FORTRAN program or read from a data file. The former method was used 
here because data files took a great amount of file space due to the large volume of 


values. 
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Figure 3.4 2-D circular plot (DISSPLA). 


Two-dimensional graphs were initially tned. By holding one parameter constant 
and with the use of some array variables, plots of directivity versus angle were 
obtained. Examples of two kinds of displays are shown in Figures 3.3 and 3.4. 
Programs are attached in Appendix B. 

To predict the real shape of a directivity pattern in three-dimensions, three- 
dimensional graphics had to be developed. In trying to get some satisfying three- 
dimensional plotting, all methods in DISSPLA were tried. They are SURFUN, and 
SURMAT for regular matrix, scattered points, and three-dimensional vector drawing. 

SURFUN is a subroutine that draws the surface from points in the z direction in 
conjunction with a grid in the x-y plane on the condition that z= f(x,y) is given. One 
constraint that was faced in applying this function was that there was more than one Z 


value for a specific x-y coordinate. 


ZZ 
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Figure 3.5 3-D plot(SURMAT for scattered points). 


The next trial was SURMAT for a regular matrix (for example, a rectangular 
three-dimensional lattice of points). For the same reason as in applying the SURFUN 
subroutine, this also turned out to be non-applicable even though all data points were 
translated into a regular matrix. 

The next option tried was SURMAT for scattered points (for example, a three- 
dimensional irregular cluster of points). As seen in Figure 3.5 it also was unsuccessful 
due to more than one point on one specific coordinate of the x-y plane. The program 
is listed in Appendix C for reference. 

The last option, the three-dimensional vector drawing method gave us a-more 
understandable display as shown in Figure 3.6. However, it did not show the beam 
pattern clearly because all the lines were visible, even from different view angles. At 


this stage the problem was the removal of hidden lines. Unfortunately, that was not 
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possible because the starting and ending points of the hidden lines could not found. 
The program is listed in Appendix D. 
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Figure 3.6 3-D plot(vector drawing method). 


Consequently, with the consideration of all the constraints mentioned above, 
another type of format was made to see the details of the beam pattern. Now, a plot 
was made using the function SURMAT about the 8-® plane instead of the x-y plane. 
This eliminated the previous multiple value problem. The results showed the most 
desirable presentation among all those tried. As shown in Figure 3.7 it gives a clearer 
appreciation of the directivity pattern in conjunction with 8 and ®. The program is 
listed in Appendix E. | 

We can relate this three-dimensional graph (Figure 3.7) to the two-dimensional 


graphs (Figures 3.3 and 3.4) for interpretation. As shown in Figure 3.3 the x-axis is 
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the angle and the y-axis is the directivity which is expressed in dB. Figure 3.4 is the 
circular graph with angle working counterclockwise in 10° increments and each circle is 
an increment of 10 dB with the center of - 40 dB. Figure 3.7 is the three-dimensional 
graph from the view point of 20° from the x-z plane, and 60° from the x-y plane. Here, 
the x-axis is @, the y-axis is 8 and the z-axis is the beam pattern (dB). From Figure 
3.3 four major points can be found at approximately 50°, 130°, 230° and 310° of 6 
with ® fixed at 126°. From Figure 3.7 the same feature can be determined at 50° of 6, 
126° of ® as seen in Figure 3.3. Since Figure 3.3 was made by holding ® constant at 
126°, one can look along the P = 126° line in Figure 3.7 and note the null at @ = 
90°, the maximum at 8 = 50°, and the relative maximum at 8 = 0°. Note the heavy 


mes In Figure 3.7. 


DIRECTIVITY PATTERN 
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Figure 3.7 3-D plot(SURMAT in 9-® plane). 
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From the circular graph (Figure 3.4) four major points can be found approximatly at 
40°, 140°, 220° and 320° of ® with @ fixed at 35°. In the three-dimensional graph one 
of these major points is seen approximately at 35° of @, 40° of M and the other major 
point is seen approximately at 35° of 8, 140° of M. Comparison is made easier when 
we see Figure 3.8 and Figure 3.9 which were drawn with the surfaces cut at D = 126° 
and @ = 35° respectively. Note the shift in axes between the two figures. All these 


points of two-dimensional graphs correspond to those of the three-dimensional graph. 


DIRECTIVITY PATTERN 
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Pe=0 P3=0 P4=0 
A=L B=1. C=L 
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Figure 3.8 3-D plot((SURMAT in 8-® plane). 


Inversely, every point of the three-dimensional graph can be read in two-dimensional 
graphs to see detailed values whenever needed. 


Finally, a trial was made to draw a constant value plot of the beam pattern. The 
first approach was to derive the functions of ® = £ (8) or ® =. f(@) and plot these on 
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DIRECTIVITY PATTERN 


K=10.5 

Di=05 D2=05 
D38=05 D4=05 
VIEW ANGLE=+20,+60 
P2=0 P38=0 P4=0 
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\ Beso f 


SS 
AEP hy 
tee, 
PATA 


ps 


SS 


HH 
wi 
yi /- 





Figure 3.9 3-D plot(SURMAT in 9-® plane). 


the 0-® plane. However, failing to derive these relationships for the general four 
element case made this method of making a contour plot of a fixed F value impossible. 
Instead, a modified version of the existing 6-® program was used. It sorted on a small 
range of values about the desired value of F to find the corresponding 8-® coordinates 
which were then plotted as in Figure 3.10. This type of plot is included throughout the 
study. 

With these forms of display and the speed of the FORTRAN version of this 
program we can now use this method to study the directional and frequency properties 


of shaded and phased simple arrays. 


Pay 


TRACK (—5 DB) 


a0. 1000 1200 ~~+*1400 
PHI(DEG) 


Figure 3.10 contour plot (1 RACK). 
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Z = —5 DBK 

K=17.568 

D1=0.3 D2=0.3 
D3=0.15 D4=0. 

P2= 0 P3=180 P4= 0 
A=1. B=1. C=0. 





IV. COMPARISON WITH SAMPLE CASES 


In making comparisons, some simple parameters were used to make the 
comparisons easier. There are three cases included in this chapter. Each case has a 
unique set of input data, and computer generated and experimental plots. In all cases, 
a two points were picked and hand calculations made to verify the computer output. 
An example of the hand calculations is given for the first case only. Comparisons were 
also made between hand or computer generated values and the two-dimensional polar 
plot, between the theoretical two-dimensional plot and the experimental data plot, 
between the two-dimensional plot and the three-dimensional plot, and finally between a 
. fixed value of the beam pattern (a “track”) and its three-dimensional plot. 

While in general, 6(®) or (8) for a constant F could not be obtained 
analytically, in a few cases, included below, it was possible to obtain the functional 


dependencies. These are discussed in the following cases. 


A. CASE 1 

When kD, = kD, ekD = KkDy = 5-2, 
~, = 9,93; = na i, 
A=B=C=1 

equation (2.37) can be simplified as 


H(9,®) = a |cos(5.275 sinOsin®) - cos(5.275 sin@cosP)}. (4.1) 


Using the above input data with values of 8 = 40° and M = 0° we get H(O,M) = 
0.984566 and F (8,M) = -0.135104 from the computer output as seen in Appendix G. 
If we use the above simplified equation and hand-computation, we get H(0,M) = 
0.98457 and F(8,M) = -0.135068 . 

In the same way, if 8 = 90° and M = 0° then we get H = 0.233302, and F = 
-12.641622 from the computer output and H = 0.23330, and F = -12.6417 from the 
hand-calculation. These hand-calculated values are almost the same the the computer 
output. Comparisons can also be made between these selected values and the two- 


dimensional plot, Figure 4.la. Additionally, the above two-dimensional plot 
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corresponds well in shape with Figure 4.1b which is from LCDR John Butler's 
experiment [Ref. 10]. ; 

If we compare the shape of the two-dimensional plot (Figure 4.la) with that of 
the three-dimensional plot, Figure 4.lc, we see that beam pattern of the two- 


dimensional graph is identical in the ® = 0° plane of the three-dimensional plot. 
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PHE 0 DEGREE 
K=10.55 

Di=0.5 D2=05 

D3=0.5 D4=0.5 

Pe=Q P3=180 P4=180 
A=1. B=L C=1. 
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Figure 4.la 2-D plot from the computer output. 


A two-dimensional plot of equation (4.1), of F values between -10.9 dB and -10.0 
dB results in Figure 4.1d. For F values less than -90.0 dB we get Figure 4. le. 
In the case where F < -90dB,H < 10°79/20 which is small enough to be assumed to 
be 0, |c,| ~ |c,] or 8 = nt n = 0,1,2---, and equation (2.40) can be tan® = 1 or -l. 


From these resulting equations we can extract ® = 45° + n7/2,n = 1,2 ---. As seen 
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in Figure 4.le, when F < -90 dB the asymptote matches the nulls in the 0-® 
relationship of Figure 4.1c. Note the nulls at @ = 0°, 180° for all ®, and M = 45°. 
These figures of constant dB level aid in the interpretation of Figure 4.lc. The 
modified P3D program (TRACE) 1s included in Appendix F. 


CO ie OE SEAM AI TERN CdB) 


R=3.20 M  kd=10.55 FREO=960 Hz 
PHASE=0, 180, 0,180 Deg AMPL= 0.265 Vac ANGLE=0 Deg 





Figure 4.1b 2-D plot from the experiment. 
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Figure 4.1c 3-D plot from the computer output. 
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Figure 4.1d Track (-10 dB). 
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Figure 4.le Track (-90 dB). 
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Be CASE 2 
In this case, a triplet array is investigated. Referring to Figure 2.3, the fourth 

element is turned off and its Opposing element is moved toward the center of the array 
So its distances is half that of the two elements. Also, the source 2 is 180° out of phase 
with the other two elements. 
When kD, = kD, = 5.2758, 

eee 2-001 7, 40, = 0, 

0, = 7,9; = 9, = 9, 

A=B=1,C=0 
the equation (2.37) can be simplified as 


I 
H(0,@) = ae {[cos(2.6379 sinOcos®)]* + [2 sin(5.2758 sinOsin®) 
- sin(2.6379 sinOcos®)]*} 4. (4.2) 


If values are calculated from the equation (4.2) as in case 1, then comparisons can be 
made between these calculated values and the two-dimensional plot as seen in Figure 
4.2a (for 8 = 0°, 20° and 90°, and D = 90°, we get F = -9.5424, -2.7431 and -3.6776 
respectively from hand-computation.). Additionally, the two-dimensional plot 
corresponds well in shape with Figure 4.2b from the experiment (Ignore the ordering of 
the relative phases in this figures; it is the consequence of a different labeling 
convention). 

If we compare the shape of a two-dimensional plot (Figure 4.2a) with that of the 
three-dimensional plot (Figure 4.2c), we see that beam pattern of the two-dimensional 
graph is identical in the ® = 90° plane of the three-dimensional plot. 

For a contour plot of this case, a two-dimensioal plot of equation (4.2) of F 
values between -5.9 dB and -5.0 dB results Figure 4.2d. For F values between -10.9 dB 
and -10.0 dB, we get Figure 4.2e. For F values less than -35.0 dB, we get Figure 4.2f- 
From equation (4.2), we get F = -35.876 for 8 = 143° and ® = 10° from hand- 
calculation. This point is seen in Figure 4.2f. Successive examination of Figure 4.2d 
through 4.2f shows how the contour lines migrate as lower and lower values of F are 
examined. These figures of constant dB level aid in the interpretation of Figure 4.2c. 
The modified P3D program (TRACE) is included in Appendix F. 
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Figure 4.2a 2-D plot from the computer output. 
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Figure 4.2b 2-D plot from the experiment. 
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TRACK (—5.0 DB) 
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Figure 4.2d Track (-5 dB). 
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Z = —5.0 DB 

K=17.586 

D1i=0.3 D2=0.3 
D3=0.15 D4=0. 
P2=180 P3=0 P4=0 
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TRACK (10.0 DB) 
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Figure 4.2e Track (-10 dB). 
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180.0 
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TRACK (—35.0 DB) 
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Figure 4.2f Track (-35 dB). 
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Z = —35.0 DB 
K=17.586 

D1=0.3 D2=0.3 
D3=0.15 D4=0. 
P2=180 P3=0 P4=0 
A=1. B=1. C=0. 





C. CASES 
In this case, a doublet array is investigated. Referring to Figure 4.2a, this ts the 

same as the first case except that the third and fourth elements are turned off and all 
elements are in phase. 
For kD, = kD, = kD, = kD, = 5.275, 

9, = ©, = Oy = Q, 

A=1,B=C#=0 
equation (2.37) can be simplified as 


H(6,M) = |cos(5.275 sin@sin®)| 


As in cases 1 and 2, for @ = 0°, 35°, 90° and D = 90°, we get F = 0, -0.058 and 
-5.459 respectively from hand-computation. These points are found in Figure 4.3a and, 
the two-dimensional plot again corresponds well in shape with Figure 4.3b from the 
experiment. 

If we compare the shape of a two-dimensional plot (Figure 4.3a) with that of the 
three-dimensional plot (Figure 4.3c below), we see that the beam pattern of the two- 
dimensional graph is identical in the ® = 0° plane of the three-dimensional plot. 

For a constant value of H(8,@), we can follow the same procedure to get the 
equation for contour plots as in case | and 2. The resultant figures are Figures 4.3d, 
4.3e and 4.3f. Figure 4.3d is drawn for the F values between -5.9 dB and -5.0 dB and 
Figure 4.3e is for the F values between -8.9 dB and -8.0 dB. Figure 4.3f is for F values 
less than -35.0 dB. 
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Figure 4.3a 2-D plot from the computer output. 
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DIPOLE BEAM PATTERN (dB) 


~ 
“ _- m0. ~~ 
“ Re N\ 
co gem 
, Tages 2 Pd ~ S,. 
oe \ \ 
is A : Btaxe “=C ~_ \ aN 
y, fo ™~ / 
i : . ee 
/ 7 o= \ \ 
| / / \ \ \vire 
=, 9 
\ , 
\ / a 
\ ‘ ‘ = y / a 
\ / /* 
| \ ar: 2 
“ \ a 
ee a a / ay, 
‘ NS _ aes : 
\ 7 ee ey 
a ee Zo tere Sf 
\itee . ~ a = ss wv / 
\ _ fa 
NX : “ 
~ os we - 


R-3.20 YH kd=10.55 FREQ=960 Hz 
PHASE=0.0 Deg AMPL=6.53 Vac ANGLE=90 Dag 


Figure 4.3b 2-D plot from the experiment. 
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Figure 4.3c 3-D plot from the computer output. 
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Figure 4.3d Track (-5 dB). 
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Figure 4.3e Track (-8 dB). 
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Figure 4.3f Track (-35 dB). 
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V. EVALUATION 


In the last chapter, we looked at three simple cases. These established the 
accuracy of the five computer programs. For example, in comparisons between hand 
calculation and computer output, both sets of values were proven nearly identical. 
Again, due to the complexity of calculations, not all values were hand-calculated. 

When a two-dimensional circular graph was compared with the selected hand- 
computed values, both groups of values were again virtually identical. As another test 
of the validitv of the five computer programs, some of the computer generated two- 
dimensional plots were compared against those from an experiment [Ref. 9: p.18]. 
Visual comparisons of the two types of plots were similar in shape, but values in dB 
were somewhat different. These differences might be caused by the normalization of 
the angular dependance of the sound pressure. Also, because of the directivity of the 
real sources, there is a slight decrease in dB level of the experimental plots when the 
array pointed awav from the receiver. However, they are small enough that the 
comparisons still seem quite reasonable. 

Next, two-dimensional versus three-dimensional graphs were compared. 
Agreement here shows the programs are working properly. By varying the range of 0 
and ® that are plotted, and also by varing the view angles, values and the shape of the 
curve can be visualized without great difficulty. 

Tracks of certain directivity values were drawn using the modified P3D program 
(TRACE) and were compared to the three-dimensional plots. The shape of the 
resulting tracks corresponded well with those of the three-dimensional patterns, again 
showing proper functioning of the programs. 

With all the comparisons made in Chapter IV, the five computer programs 
developed in this thesis were proven effective in predicting the directivity pattern in all 
directions of space with or without differences in the strength of ‘the sources, inter- 


source distances and phases. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 
The following conclusions are offered. 
e The P2DICIR, P2DPCIR, P2DIFIX, P2DPFIX programs successfully 
predict the directivity beam pattern in two dimensions. 
¢ The P3D program predicts and visualizes the directivity beam pattern in the 
three-dimension space. 
¢ They allow the presentation of complete beam patterns, which contain 


considerable information, in a succinct, easy to comprehend and interpret form. 


B. RECOMMENDATIONS 
The following recommendations and areas of further investigation is suggested. 
e Further studies for computer programs to visualize the beam pattern in the 
X-V-Z space are needed to eliminate hidden lines in three-dimensional plotting. 
e Transcendental equations need to be solved to get the contour plot for fixed 


F values. 
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BASIC PROGRAMS (HP COMPUTER) 








rs % 
*% THE OBJECTIVE OF THESE PROGRAMS IS TO CALCULATE AND * 
*% DRAW THE DIRECTIVITY PATTERN IN TWO-DIMENSION. % 
% THIS IS WRITTEN IN BASIC LANGUAGE AND PLOTTER IS HP7090A. *% 
* : % 
| WRBBRRERHHHRHEHHHEEHHERH VARIABLE DEFINITION  24362636363636969636363696065696066 20363636269630063¢ 
% * 
3% Zz : DIRECTIVITY IN cB % 
*% P2 : PHASE (DEGREE) OF SOURCE A 3% 
*% P3 ;: PHASE (DEGREE) OF SOURCE B % 
*% PG ; PHASE (DEGREE) OF SOURCE C * 
% PS : VARIABLE TO CHANGE DEGREE TO RADIAN % 
*% P6 : RADIAL VALUE OF Pl % 
% P7 RADIAL VALUE OF P2 3% 
% P8 RADIAL VALUE OF P3 % 
*% ee “WAVE NUMBER % 
* D ; INTER=SOURCE DISTANCE x 
* P ;: PHI IN DEGREE 2% 
% I :; THETA IN DEGREE * 
% PI : 3.1415917 % 
& 3% 


IIE BEBBE BBB BB HERE BEBE BRR RBBB EBB HUBBER B BEBE BBE HERE BEB HERB HBE88Ht 
HP PROGRAM FOR THE RECTANGULAR GRAPH 
IEB BBB EH EBE BEB HEE EBB BBB HRB BEEBE BBB HBR EB NBR BBB 


INPUT PARAMETERS 
5 OPTION BASE 1 
6 RAD 
10 DIM 2(360) 
20 DISP "WHAT IS THE VALUE OF KD? KD= “! KD=K*D 
30 INPUT KD 
60 DISP “WHAT IS THE VALUE OF P2(DEG)? P2= "! P2=PHASE OF A 
61 INPUT P2 
62 DISP “WHAT IS THE VALUE OF P3(DEG)? P3=" ! P3=PHASE OF B 
63 INPUT P3 
64 DISP “WHAT IS THE VALUE OF P3(DEG)? P4G= " ! PG=PHASE OF C 
65 INPUT P4 
66 DISP “WHAT IS THE VALUE OF A? A= " 
67 INPUT A 
68 DISP_ "WHAT IS THE VALUE OF B? B= " 
69 INPUT B 
70 DISP "WHAT IS THE VALUE OF C? C= " 
71 INPUT C 
72 DISP_ "IS P FIXED? (YES$AN=1 , NO$SAN=0) AN= 
73 INPUT AN 
74 IF AN=1 THEN DISP "WHAT IS THE VALUE OF P(DEG)7 P= "ELSE GOTO 76 
75 INPUT P 
76 DISP_ "IS I FIXED? (YES$ANS=1 , NOSANS=0) ANS= “ 
77 INPUT ANS 
78 IF ANS=1 THEN DISP “WHAT IS THE VALUE OF I(DEG)? I= “ELSE GOTO 82 
NPUT I 
BZ DISP 8993636363636 36363636566 IESE IE IESE IEEE IE IE SESE IE IRIE NE IE IEIE IEEE IEEE IEEE IE IEEE IEEE 36 
83 IF AN=0 AND ANS=0 THEN GOTO 86 
84 IF AN=0 AND ANS=1 THEN GOTO 90 
85 IF AN=1 AND ANS=0 THEN GOTO 93 
86 DISP Oh 5g08 3"KxeD=" 3KD > oe 608 > 00 seen 
87 DISP "*PHASE OF A(DEG)="5P2;"PHASE OF B(DEG)="5P35"PHASE OF C{DEG)="5P4, "x" 
88 DISP WO xe0 "A=" 3A, 16 B=" 3B,"C=" 3C , "se" a GOTO 99 
90 DISP 06 sete "KBs" 5 KD, ta If DEG }="" 3I a “8 se 
91 DISP “XPHASE OF A(DEG)="sP2;"PHASE OF B(DEG)="sP35"PHASE OF C(DEG)="5P4,"x" 
92 DISP "x"5"A="5A," B="5B,"C="5C,""" 3 GOTO 103 
93 DISP 00 ygte 3;"Kx¥D=' 3$KD > ve Ie DEG att $i > te 8646 > Ob sgt 


aH 


94 DISP “xPHASE OF A( DEG J="3P23;"PHASE OF BI DEG)=""3P33;"PHASE OF C( DEG )=''3PG, te" 


95 OISP eeggre se Ass A,” B="3B,"Cs"3C, "x" a GOTO 107 

99 DISP [RRB ERE HEEB BE BH BR HERBEBEHERENEEEB RBH REE HEHER HEHE 
100 orIsP " 

101 OISP " "3" THETA" PHI" ," DIRECTIVITY" 

102 FOR P=1 TO 360 9 GOTO 110 

195 DISP (RB ROBB RRB BOB RHR RB HER HB RE IE HERO BRB RE ROE RHEE HRBEIEE 
104 DISP * " 

105 OISP “ "3" THETA*3® PHI"'," DIRECTIVITY" 

106 FOR P=1 TO 360 a2 GOTO 115 

LO? JISP 8 3R6RB SBR RB BOR RR ERR HEBER ROBE HHO BRB IEEE EERE RBHRE! 
108 OISsP “ " 

109 DISP " "3" THETA™;" PHI" ,"* OIRECTIVITY” 


aaa er ee eee awe ae SSP SPF 2 SSP SP FF SPF P@Ws* SF as Bee Bese ase eae sae eee SOF BFF SVWs a @aasaaaaesaaa 2 waa 


110 FOR I=1 TO 360 

115 P5=PI /180 

116 P6=P2*P5 2 P7=P3¥PS5 Dd PB=PGxPS 

117 Q=P-90 

118 Q1=SIN (QxP5} 

119 I1=SIN (I*P5) 

120 Q2=SIN (P*P5) 

121 W1=KOxI1*Q2/2 

121 W2=KOxI1*Q1/2 

127 W=COS (W1)4AxXCOS (WN14P6)4+B*¥COS (WN2-P7}4CxCOS (WN24+P8) 
128 G=SIN (W1)-A¥SIN (W14P6)4B*¥SIN (WN2-P7}-C¥SIN (N24P8) 
130 H=SQR (W¥¥24G6%*2 )/(14A4B4C } 

140 Z(P)=20*LGT (H) 

143 T=20*LGT (ABS (SIN (.5%*KOx*I1))}) 

150 DISP USING 153 3 I,P,Z(P) 

153 IMAGE 5xX,30.0,5X,30.0,15xX,S000.000,15xX,S000.000 
' IF 2Z(1)<-20 THEN I=I+.1 3 GOTO 119 

156 ! I= IP(I} 

157 ! IF Z(P)<-20 THEN P=P+.1 9 GOTO 117 

158 ! P= IP(P)} 

160 IF AN=1 ANDO ANS=0 THEN GOTO 163 ELSE GOTO 161 
161 IF AN=0 AND ANS=1 THEN GOTO 164 ELSE GOTO 163 
162 NEXT ee 96 


SE[Eee GS ee eee 2 2 S22 2 SP 2 O28 SFP SPSS @FISs* FB SFISFQ_H @ Fs BI FB FFI PWSs*Q eo BPs BF Bers FT@W@SsISP@ee2 2 ea eg as ese 8 @ @ = 


170 PLOTTER IS 705 

180 CLEAR 

190 GCLEAR 

200 GRAME 

210 LOCATE 15,145,10,80 

220 SCALE 0,390,-240,0 

230 CLIP 0,360,0,-240 

240 AXES 5,55,05-240,10,1052 
250 LINE TYPE 1 

260 GRIO 20,-20 

270 FRAME 

280 LORG 6 

290 CSIZE 4 

300 FOR I=0 TO 360 STEP 40 
310 MOVE I,-241 

320 LABEL I 

330 NEXT I 

340 LORG 8 

350 FOR I=-240 TO O STEP 40 
360 MOVE -1,I 

370 LABEL I 

380 NEXT I 

590) CSIZE TS 

400 LORG 5 

410 MOVE 180,40 

420 CSIZE x 

430 PEN 2 

440 LABEL “PHI VS OIRECTIVITY" 
450 LORG 5 

460 MOVE 180,20 

470 CSIZE 4 

480 LABEL “(THETA = 45 DEGREE)" 
490 CSIZE 5 

500 LORG 5 

510 MOVE 180,-26- 

520 CSIZE 4 


a2 


530 


540 


550 


570 


LOIR 90 


LABEL "PHI( DEG)" 


MOVE -30,-120 
560 DEG 


LDIR 90 


580 CSIZE 4 


590 


600 
610 
620 
630 
640 
650 
700 
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CSIZE 5 
MOVE 0,0 


FOR P=1 TO 360 
DRAW P>Z2(P) 


NEXT P 
PEN 0 
END 
7392 


LABEL “DIRECTIVITY( dB)" 


IBB EUR B08 8888888888888 8881888 BBE BREE IEE BBE BIRR IRI III 
HP PROGRAM FOR CIRCULAR GRAPH 
JOB EBB EBB BBR BBR IEETEEHEBESESIEDIE IESE IEEE IEE 


2&2Geee@ee2eeee eee eeee ee Se Se SO Se SSF SF SF SF SF FF FF SFP SF FB SF SF SF SFT SF Oe OF SF SF SF FB SF FSF GVWea Sere ee Oe was &S = @& & 


5 OPTION BASE 1 
6 RAD 


10 DIM X( 360) 

ll DIM Y(360) 

12 DIM 2(360) 

20 DISP "WHAT IS THE VALUE OF KD? KD= "“ !? KD=K*D 

30 INPUT KD 

60 DISP “WHAT IS THE VALUE OF P2(DEG)? P2= “t P2=PHASE OF A 
61 INPUT P2 

62 DISP “WHAT IS THE VALUE OF P3(DEG)? P3= " ! P3=PHASE OF B 
63 INPUT P3 

64 OISP "WHAT IS THE VALUE OF PG(DEG)? Pa= " ! PG=PHASE OF C 
65 INPUT PG 

66 DISP “WHAT IS THE VALUE OF A? A= "™ 

67 INPUT A 

68 DISP "WHAT IS THE VALUE OF B? B= "™ 

69 INPUT B 

70 DISP “WHAT IS THE VALUE OF C? Cs " 

71 INPUT C 

72 DISP "IS P FIXED? (YES;AN=l1. » NOsAN=Q) AN= *“* 

73 INPUT AN 

74 IF AN=1 THEN DISP “WHAT IS THE VALUE OF PI(DEG)? P= “ELSE GOTO 76 
75 INPUT P 

76 DISP “IS I FIXED? (YES;ANS=1 , NO;ANS=0)} ANS= “ 

77 INPUT ANS 

78 IF ANS=1 THEN DISP “WHAT IS THE VALUE OF I(DEG)}? I= “ELSE GOTO 82 
79 INPUT I 

B2 DISP 19363635636 26 2 IEE REI SHI IER IE ELSIE ER RT OE 


83 IF AN=0 AND ANS=0 THEN GOTO 86 
84 IF AN=0 ANO ANS=1 THEN GOTO 90 
85 IF AN=1 AND ANS=0 THEN GOTO 93 


86 DISP Hes KOH" SKN," ae, set 

87 DISP "xPHASE OF A(DEG)="3P23"PHASE 
88 DISP Toe ae eae B= to5 G=" 3G, "3" 
90 DISP "*"3"KxO="5KD,"  I(DEG)="51," 
91 OISP "*PHASE OF A(DEG)="5P25"PHASE 
92 DISP etna Sg B="3B,"CH"3C, ee" 
93 DISP "x"5"KxO="5KD," I(OEGI="5I 5" 
94 DISP "*PHASE OF A(DEG)="5P25"PHASE 
95 DISP Cet Re tA 5 B="3B,"Ca"3C, x" 
99 DISP 

100 DISP " " 

101 DISP “ "3" THETA"s" | PHI," 
102 FOR P=1 TO 360 9 GOTO 110 

103 DISP 

104 DISP " " 

MeSepISPo* Ss" THETA"; PHI") 
106 FOR P=1 TO 360 2 GOTO 115 

107 DISP "ieee 

108 DISP " " 

109 DISP " "3" THETA'3" = PHI"," 


OF BC DEG J="3P33;"PHASE OF C(I DEG )="3PG, *x* 
2 GOTO 99 
68 Ce zeas 
OF B( DEG J="3P33;"PHASE OF C(DEG }="3PG, "x" 
2 GOTO 103 
oe > Oh yee 
OF BIDEG)="3;P33"PHASE OF C(IDEG }="3P4, "x" 
2 GOTO 107. 

FIORE" 


DIRECTIVITY" 


“RBI REE GREE K BIG OSH GEESE bi IE KI IS PE 


DIRECTIVITY" 


DIRECTIVITY” 


110 FOR I=1 TO 360 
115 P5=PI /180 


S35 


116 P6=P2*P5 3 P7=P3xP5 9 P&=P4xP5 
117 Q=P-90 
118 Q1=SIN (Q*P5) 
119 I1=SIN (I¥P5) 
120 Q2=SIN (P¥P5) 
121 W1=KD*¥I1*Q2/2 
121 W2=KD¥I1*Q1/2 
127 W=COS (W1)+AxCOS (W1+P6 34+BxXCOS (W2-P7)+CxCOS (W2+P8) 
128 G=SIN (W1)-AX¥SIN (W1+P6 )+BXSIN (W2-P7)J-CXSIN (W24+P8) 
130 H=SQR (W**24+G%*2 )/(14+A+B4C ) 
140 Z(IJ=20*LGT (H)+240 
141 X(I)J=Z(IJXCOS (I¥PS) 
142 Y(IJ=Z(I *XSIN (I*¥P5) 
143 T=20*LGT (ABS (SIN (.5x*KD¥I1))) 
150 DISP USING 153 3 I,P,2(P) 
153 IMAGE 5X,SD00.000,15X,S000.0D0 
t IF 2Z(I)<-20 THEN I=I+.1 9 GOTO 119 
156! f= stPiT) 
157 ! IF 2(P3<-20 THEN P=P+.1 9 GOTO 117 
158 ! P= IP(P)) 
160 IF AN=1 AND ANS=0 THEN GOTO 163 ELSE GOTO 161 
161 IF AN=0 AND ANS=1 THEN GOTO 164 ELSE GOTO 163 
162 NEXT P 9 GOTO 96 
163 NEXT I 
164 ! NEXT P 


eS @eeeeeeeeeee eee @eee22 22 8 Ss ee Fee 82 2 2 BF SF FF SF VF SOS 288 SP SF FS ae Be 2 Fae eee 2 SO SOF SS SF SFB SF Fe SF 


170 PLOTTER IS 705 

180 CLEAR 

190 GCLEAR 

200 GRAME 

210 LOCATE 25,120,5,80 

220 SCALE -260,320,-220 »240 
230 CLIP -240,240,-240,240 
235 XAXIS 0,20,-240,240 

240 YAXIS 0,20,-240,240 

250 LINE TYPE 1 


270 FRAME 

280 LORG 6 

290 CSIZE 4 

300 FOR I=0 TO 360 STEP 40 
310 MOVE 1I;-5 

320 LABEL I-240 

321 MOVE -1,I1 

322 LABEL I-240 

330 NEXT I 

340 LORG 6 


341 MOVE 0,-5 
342 LABEL -240 
350 FOR I=-240 TO -60 STEP 60 


360 MOVE -1,I1 
370 LABEL -240-I 
So? MOVE I,-5 
372 LABEL -240-I 
375 NEXT I 


380 MOVE 0,0 

381 FOR I=0 TO 240 STEP 60 

382 FOR J=0 TO 360 

383 X=IxXCOS (J*PI /180) 
384 Y=IXSIN (J*PI /180) 
385 DRAW X>Y 

386 NEXT J 


410 MOVE 0,310 

420 CSIZE 8 

430 PEN 2 

440 LABEL "THETA VS DIRECTIVITY" 
450 LORG 50 

460 MOVE 0,285 

‘470 CSIZE 4 

480 LABEL "(PHI = 90 DEGREE)" 
490 CSIZE 5 

500 LORG 5 

510 MOVE 300,0 

520 CSIZE 4 

530 LABEL "O(DEG)" 

550 MOVE -300,0 

580 CSIZE 4 

590 LABEL "180( DEG)" 


54 


610 MOVE 0,0 

620 FOR I-1 TO 360 

630 DRAW X(I),Y(I) 
640 NEXT I 

650 PEN 0 

700 END 

111174 


a5 


APPENDIX B 
P2DIFIX, P2DPFIX, P2DICIR AND P2DPCIR 





: 
3 





* * 
% THE OBJECTIVE OF THESE PROGRAMS IS TO CALCULATE AND % 
* DRAW THE DIRECTIVITY PATTERN IN TWO-DIMENSION. % 
% THE DRAWING METHOD IS THE SUBROUTINE CURVE OF DISSPLA. % 
* I(U) AND P(S!) CAN BE USED INTERCHANGEABLY FOR * 
* P2DPFIX AND P2DICIR PROGRAMS. % 
x * 
VARIABLE DEFINITION 9 206888BRBBBBRBEHREBREHRERHEHE 
* * 
* cn: DIRECTIVITY IN cB x 
* Pl PHASE (DEGREE) OF SOURCE A a 
* P2 : PHASE (DEGREE) OF SOURCE B % 
% P3 : PHASE (DEGREE) OF SOURCE C % 
* PG : VARIABLE TO CHANGE DEGREE TO RADIAN * 
* P5 : RADIAL VALUE OF Pl % 
% P6é RADIAL VALUE OF P2 x 
* P7 RADIAL VALUE OF P3 * 
% K : “WAVE NUMBER * 
% Dl : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE 1 * 
* D2 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE A % 
% D3 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE B * 
% DG ;: DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE C % 
% P ; PHI IN DEGREE * 
% I : THETA IN DEGREE % 
* PH : PHI IN RADIAN x 
* TH : THETA IN RADIAN % 
x PHI : DUMMY VARIABLE OF PHI * 
% THETA : DUMMY VARIABLE OF THETA % 
% PI : 3.14¢15917 * 
*% e * 
JBRRREHRRRRHRBAHRHHHBHHE VARIABLE DECLARATION 280RRRRRHRHRHHRHRRHHHHHHEHE 
Cc 
C 


9993303 EIEIO EE FEE IEEE FE ESE DE SESE FE FEFE SESE SE FEE SESE FEE FE SEE FE FE FEE EE TE 


PROGRAM FOR THE TWO-DIMENSIONAL RECTANGULAR GRAPH 
PBB RHEUN IBA OG BOBBIE IBF 


Cc 
DIMENSION PHI(361}),THETA(361),2Z( 360} 
REAL P4,P5,P6>,P7,W1,W2,WN3 »W4,N»G>H»K»D1,D2,03,D4¢,PH,TH»PI>,a,bs¢ 
: INTEGER P,P1,P2,P3,1 
ee 
; INPUT PARAMETERS 
Cc 
PI=3.1415917 
K=5.5 
D1=1. 
D2=1. 
D3=1. 
DG=1. 
Pl= 180 
P2= 0 
P3= 180 
A=l. 
B=l. 
c=l. 
5 I = 25: USE P FOR P2DPFIX 
Cee m meee e eee w woe w we ewes coe ewww ween$ eww eeeseeecoeeeeooe ow ee oooeocoecoocecns 
e MAIN PROGRAM 
Cc 
P4=PI/180 
P5=P1¥P4 
P6=P2*P4 
P7=P3¥P4 
¢ DO 100 I=1,360: MOVE C TO LINE NO. 76 FOR P2DPFIX 


DO 200 P=1,360 


56 


TH=I¥#P4 
PH=P¥P4 
W1=K¥D1*SIN( PH J¥SIN(TH ) 
W2=K¥D2*SIN( PH )¥SIN( TH )_P5S 
W3=K¥D3*SIN( TH }*COS( PH }-P6 
W4=K*¥D4*SIN( TH J®COS( PH )+P7 
W=COS(W1 3 +A*COS( W2 3 +B¥COS( WS 3+C¥COS( WG } 
G=SINCUW1 J-AXSIN(OW2 )-BXSINC WS J4CXSIN( WG J 
H=SQRT (WHN+G#G )/( 1+A+B4C J 
Z(P )=20*ALOG10(H} 
IF(Z(P).LE.-40.) THEN 

Z(P )=-40 


SE 
Z(P)=Z(P } 
ENO IF 
THETA(P )=I 
PHI(P )=P 
WRITE(8,1) PHI(P),Z( 
FORMAT( 2X%,2(E10.3,5X 


— TW 
—w 


200 CONTINUE 
100 CONTINUE 


ME DBUF 

COMPRS 

TEK618: MOVE C TO ORAW USING TEK618 
PHYSOR(1. 5,1. } 

PAGE(15.5,12.) 

AREA2D(10.5,9.1) 

FRAME 

COMPLX 

HEIGHT( .2) 

XNAME( 'PHI( DEG) $',10) 

YNAME( 'DIRECTIVITY(0B8)$',15) 
YTICKS(8 ) 

XTICKS(8 ) 

GRAF(0,40,360,-40,4%,0) 

GRID(2,2) 

CURVE (PHI ,Z,359,0) 

MESSAG( ‘THETA=25 DEGREE',15,11.,8.5) 
MESSAG( ‘'K=5.5',5,11. »8.2) 
MESSAG('O1=1. O2=1.',11,11.,7.9) 
MESSAG( 'D3=1. DG=1.',11,11.5,7.6) 
MESSAG( 'P2=180 P3=0 P4=180',18,11.,7.3)3 
MESSAG( 'A=1. B=l. C=1.',14,11.5,7.0) 
RESET( ‘HEIGHT ' } 

HEIGHT( .4) 

MESSAG( ‘PHI VS. DIRECTIVITY' ,19,3.1,9.8) 
RESET( 'ALL' } 

NOCHEK 

ENDGR(0) 

PHYSOR(11.85,7.6) 

AREA2D(3.5,2.5) 

FRAME 

ENDPL(0 } 


PROGRAM FOR THE CIRCULAR TWO-DIMENSIONAL GRAPH 
JOBE BRB REE BERBER BOGE RE 


FHA IONE 


% 


OIMENSION PHI(360),R(360 ),THETA( 360) »>2(360),X(4),Y(4) 


REAL P4,P5,P6 »P7»W1 NZ WS NG W>GH»K»D1,02,D3,D4,PH,TH»PI ,asb5c 
INTEGER P,P1,P2,P3,I,Q 
Comm m www eww ewww www wwe wwe cow e ewww wee on cee we eww www owt eee we wee e owe eee ooo wee 
: INPUT PARAMETERS 
C 

PT=3.1415917 

K=5.5 

D1l=1. 

D2=1. 

D3=1. 

DG=1. 

Pl= 180 

P2= 0 

P3= 180 

A = ce 


a7 


4 


P4=PI/180 
P5=P1*xPq 
P6=PZ2*PG 
P7=P3*PG 
DO 100 I=0,;360: MOVE C TO LINE NO. 174 FOR P2DICIR 
DO 200 P=0,360 
TH=I¥PG 
PH=PxPG 
W1=K*¥D1*SIN( PH )¥COS( TH ) 
W2=K¥DZ2¥SIN( PH )*XCOS( TH )+P5 
W3=K*D3*SIN( TH )¥COS( PH )-P6 
WG=K¥DGXSIN( TH )¥COS( PH )4+P7 
W=COS(W1 )+A¥COS(W2 )+B*COS( W3 )+C¥COS( WG) 
G=SIN(W1 J-AXSIN(WN2 )+BXSIN( W3 )-CXSIN( WG ) 
H=SQRT (W¥N4+G¥G )/(14A+B4C ) 
Z(P )=20%ALOG10(H ) 
IF(Z(P).LE.-40. ) THEN 
2(P)=0. 
ELSE 
Z2(P )=40.4+2Z(P ) 
END IF 
THETA(P )=I 
PHI(P)=P 
WRITE(8,1) PHI(P),2Z(P) 
FORMAT( 2X,2(E10.3,5X)) 
200 CONTINUE 
100 CONTINUE 


CALL MEDBUF 

CALL COMPRS 

CALL TEK618: MOVE C TO THE UPPER LINE WHEN DRAWING WITH TEK618 
CALL PHYSOR(1.5,1.) 

CALL PAGE(15.5,12.) 

CALL AREA2D(2. >8. ) 

CALL COMPLX 

CALL HEIGHT( .2) 

CALL POLAR(3.1415917/180. »10. »%. 5%. ) 

CALL GRID(1,1) 

CALL CURVE(PHI,Z,360,0) 

CALL MESSAG( 'THETA=25 DEGREE',15,9.7>,7.7) 
CALL MESSAG( "'K=5.5' 959.797.-4%) 

CALL MESSAG( ‘D1=1. D2=1.°,11,9.7,7.1) 

CALL MESSAG( 'D3=1. 04=1.°,11,9.7,6.8) 

CALL MESSAG( 'P2=180 P3=0 P4=180' ,18,9.7,6.5) 
CALL MESSAG('A=1. Bl. C=1.'514,9.7,6.2) 
CALL MESSAG(‘'0',1,3.9,8.1) 


CALL MESSAG( '-10',3,3.7,7.1) 
CALL MESSAG( '-20' ,3,3.7,6.1) 
CALL MESSAG( '-30',3,3.7,5.1) 
CALL MESSAG( '-40' ,3,3.7,54.1) 
CALL MESSAG( ‘O(DEG)'»,6,7.65,3.9) 
CALL MESSAG( '180°,3,0.,3.9) 


CALL RESET( "HEIGHT ') 

CALL HEIGHT(.4) 

CALL MESSAG('PHI VS. DIRECTIVITY',»,19,3.0,9.5) 
CALL RESET( 'ALL*) 

CALL NOCHEK 

CALL ENDGR(O) 

CALL PHYSOR(10.95,6.8) 

CALL AREA2D(3.5,2.5) 

CALL FRAME 

CALL DONEPL 


58 


APPENDIX C 
PROGRAM FOR THE 3-D GRAPH USING SURMAT ABOUT X-Y PLANE 








INFORMATION  36363636389090369636900696309036563€ 





a * 
% THE OBJECTIVE OF THIS PROGRAM IS TO CALCULATE AND DRAW % 
* THE DIRECTIVITY PATTERN IN THREE-DIMENSION. % 
% THE DRAWING METHOD IS THE SUBROUTINE SURMAT OF DISSPLA. ¥ 
* * 
FBO BREeHeHeeHRReeH «VARIABLE DEFINITION 23393RB0B0HHHHHHRHHHERHRHIIE 
* % 
% Z : DIRECTIVITY IN cB * 
* 2MM ;: DOUBLE MATRIX TO GET X-Y PLANE * 
% 2M : DUMMY VARIABLE OF 2 * 
* Pl : PHASE (DEGREE) OF SOURCE A % 
* P2 : PHASE (DEGREE) OF SOURCE 8 * 
* PZ : PHASE (DEGREE) OF SOURCE C * 
* PG : VARIABLE TO CHANGE DEGREE TO RADIAN % 
* P5 : RADIAL VALUE OF Pl * 
% P6 : RADIAL VALUE OF P2 % 
* P7 : RADIAL VALUE OF P3 * 
* K : WAVE NUMBER * 
* Dl : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE 1 * 
* D2 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE A % 
* D3 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE B * 
% D4 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE C * 
¥ PHI : PHI IN DEGREE * 
* I : THETA IN DEGREE * 
* PH : PHI IN RADIAN * 
* TH : THETA IN RADIAN % 
* PI : 3.1415917 % 
*% * 


FOSRHHHREROHBHHHEHHHHeHE VARTABLE DECLARATION 2228 RHRHRHRHHRRORHE 
DIMENSION 2(8100),ZMM( 90,90) 
REAL P»>TsR»XsY»Z2Msasboc 
REAL PG,P5,P6,P7,W1,W2Z,WN3 »WG4,WN,CG,H»K»D1,02,03,04,PH,TH,PI 
INTEGER PHI,P1,P2,P3,I 


CALL COMPRS 

CALL TEK618: MOVE C TO THE UPPER LINE WHEN DRAWING WITH TEK618 
CALL RESET (3HALL) 

CALL PAGE(15.,9. ) 

CALL PHYSOR(1.,3.) 

CALL SCMPLX 

CALL AREA2D(14.,8. ) 

CALL X3NAME(1H ,1) 

CALL Y3NAME(1H 51) 

CALL Z3NAME(1H ,»,1) 

CALL VOLM3D (17. 517.517. ) 

CALL VUANGL(+4+20.,+420. ,500. ) 

CALL INTAXS 

CALL ZAXANG (-90. ) 

CALL GRAF3D(0.,10.,80.,0.,10. ,80.,0.,10. ,40. ) 
CALL BLSUR 

CALL BGNMAT( 90,90) 


PI=3.1415917 
K=PTI 
D1l=1. 


ao 


Dp3=1. 
D¢=1. 
Pl= 180 
P2= 0 
P3= 180 
A=il. 
B= 0. 
Cc = 0. 
c 
Ca a eee 
c MAIN PROGRAM 
 é QE SS ee eee e sree SF eee eneeeenvowvwwsowc wesw sooo Ss eee e cance esas anwowwwwww ooo seeeeeoooe 
c 
PG=PI/180 
P5=P1*PG 
P6=P2*PG 
P7=P3*PG 
DO 200 I=-90,0 
TH=I¥PG 
DO 100 PHI=1,360,% 
PH=PHI*PG 
W1=K¥D1L¥SIN( PH J¥SIN(TH ) 
W2=K¥D2*SIN( PH J¥SIN( TH J-P5 
W3=K¥D3¥*SIN( TH )*COS( PH J-P6 
WG=K¥DGXSIN( TH )*COS( PH )+P7 
W=COS(W1 }+A*¥COS(W2 }+BxCOS( W3 )+CxXCOS( WG ) 
G=SIN(W1 ) ~A¥XSIN(W2 ) -BXSIN( WZ ) ¢C¥SINCWG ) 
H=SQRT ( W*#H+GG )/(14+A4+B4C ) 
Z( PHI J)=20*ALOG10(H } 
IF(Z( PHI }.LE.-40.)} THEN 
Z( PHI )=0. 
ELSE 
Z( PHT )=Z( PHT ) 
END IF 
R=40.+2Z( PHI ) 
P= FLOAT(PHI ) 
T=FLOAT(TI ) 
X=R*¥COS( TH )¥COS( PH )+40. 
Y=R¥COS (TH )¥SIN( PH )+40. 
=-R¥SIN( TH ) 
WRITE(8S,1) Xs¥,2ZM 
1 FORMAT( 2% ,3(E10.3,5xX)) 
CALL GETMAT(X,Y,Z2M5,1;,0) 
100 CONTINUE 
200 CONTINUE 
. 
: A PROGRAM FOR PLOTTING BY SHERPAOR BY TEK618 
C 
CALL ENDMAT(2ZMM;0 ) 
CALL SURMAT (2MM,1,90,1,90,0) 
CALL DASH 
CALL MARKER(17) 
CALL COMPLX 
CALL HEIGHT(.2} 
CALL SURVIS( ‘BOTTOM’ ) 
CALL MESSAG( 'DIRECTIVITY PATTERNS',19,2.0,-.5) 
CALL RESET( 'COMPLX' } 
CALL RESET( ‘HEIGHT’ ) 
Cc 
Cc 
PROGRAM FOR THE LOWER HEMISPHERE 
Cc 
Cc 
- CALL BGNMAT( 90,90) 
: INPUT PARAMETERS 
C 
PI=3.14615917 
K=PI 
Di=1. 
D2=1. 
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D3=1. 


PG4=PI/180 
P5=P1¥*¥PG 
P6=P2*PG 
P7=P3*PG 
DO 250 I=-180,-90 
TH=I¥PG 
DO 150 PHI=1,360 54 
PH=PHI*PG 
W1=K¥D1*SIN( PH )¥COS( TH ) 
W2=K¥DZ2*SIN( PH J*COS( TH )+P5 
W3=KXD3*SIN( TH )*XCOS( PH J-P6 
NG=K*xDGXSIN( TH J*COS( PH J4+P7 
W=COS(W1 )+A¥COS(W2 )+BxCOS( WS )4CxXCOS( NG ) 
G=SIN(W1 )-AXSINOW2 J) 4BXSIN(W3 J-C¥SINOWG ) 
H=SQRT (WeWN4+GXG )/(14A4B4C J 
Z( PHI J=20*ALOG10(H ) 
IF(2Z(PHI).LE.-40.) THEN 
Z( PHI )=0. 
ELSE 
Z( PHI J=Z( PHI ) 
END IF 
R=40.+Z( PHI ) 
P= FLOAT(PHI ) 
T=FLOAT(T ) 
X=R¥COS( TH )*COS( PH )440. 
Y=RxCOS( TH )¥SIN( PH J4+40. 
ZM=R¥SIN( TH J 
WRITE(8,1) X,¥,2M 
1 FORMAT (2% ,3(E10.3,5X)) 
CALL GETMAT(X>,Y,2M,1,0) 
150 CONTINUE 
250 CONTINUE 
CALL ENDMAT(ZMM,0) 
CALL SURMAT (2ZMM,1,90,1,90,0) 
CALL DASH 
CALL MARKER(17 ) 
CALL COMPLX 
CALL HEIGHT( .2) 
CALL SURVZS( 'BOTTOM' ) 
CALL ENDPL(0) 
CALL DONEPL 
CHECKSUM=012362 
STOP 
END 
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APPENDIX D 
PROGRAM FOR THE 3-DIMENSIONAL GRAPH USING VECTOR DRAWING MET 


SO RII SUBURB INFORMATION 268299 B 8B BBE BHU OHHHGE 
x x 
*% THE OBJECTIVE OF THIS PROGRAM IS TO CALCULATE AND DRAW x 
* THE DIRECTIVITY PATTERN IN THREE-DIMENSION. * 
* THE DRAWING METHOD IS THE VECTOR DRAWING METHOD OF DISSPLA. * 
x » 
JOIBEBEBEHRBRBEBHBERRHHE «VARIABLE DEFINITION 33338020 BBBHRURRRUHRROBHHEE 
x ce 
* Za: DIRECTIVITY IN dB 3 
* XT "TO" POINT OF X IN DRAWING VECTOR % 
* YT : “TO POINT OF Y IN DRAWING VECTOR * 
3 Zr : “TO POINT OF Z IN DRAWING VECTOR 36 
* XF : “FROM POINT OF X IN DRAWING VECTOR * 
x YF : “FROM POINT OF Y IN DRAWING VECTOR % 
ue ZF : “FROM POINT OF Z IN DRAWING VECTOR % 
x Pl : PHASE (DEGREE) OF SOURCE A x 
* P2 : PHASE (DEGREE) OF SOURCE B * 
* P3 : PHASE (DEGREE) OF SOURCE C * 
% PG : VARIABLE TO CHANGE DEGREE TO RADIAN * 
* PS : RADIAL VALUE OF Pl * 
% P6 RADIAL VALUE. OF P2 2 
% P7 RADIAL VALUE OF P3 % 
% K : "WAVE NUMBER % 
* Dl : DISTANCE: FROM THE GEOMETRICAL CENTER TO THE SOURCE 1 % 
% D2 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE A *% 
% D3 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE 8B % 
*% DG : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE C % 
% PHI : PHI IN DEGREE * 
% I : THETA IN DEGREE % 
% PH ;: PHI IN RADIAN * 
* TH : THETA IN RADIAN * 
% PI : 3.1415917 * 
* * 
JABBER HHHRBHRHHHERHHE «VARIABLE DECLARATION 28208 RHRRRRRRHRHRHHERHHHEHHEIE 


REAL R5Z oXT o¥T »XF »YF »r2F oZI oa »>boc 
REAL PG,P5,P65P7,W1 WZ »N3»NG,N,G>H»K»,D1,D02,03,04¢,PH,TH,PI 
INTEGER PHI,P1,P2,P3,I 

Cc 


CIB BE BEBE BBR EBB BEEBE NBER B88 R868 RRB BBB BERBER 
DIMENSION 2Z(0:2500 )»ZMM(46>4%6 ),V( 2000) ,U( 2000) 
REAL P»T»RoXsr¥»ZM»oZZM oA oBoCoV>U 
REAL PG,P5»>P6 »P7,W1,W2 »W3 »NG,N»G»H»K»D1,D02,03 »D4,PH,TH,PI 
INTEGER PHI,P1,P2,P3,15S 


C 
ee Se Oe eS Se SS eS Se ee 
INPUT PARAMETERS 
Cc 

PI=3.1415917 

K=10.55 

D1=0.5 

D2=0.5 

D3=0.5 

D4=0.5 

P1l= 0 

P2=180 

P3=180 

A=l. 

B= 1. 

c=l. 
C 
nn ee oe ee eae oa 
s MAIN PROGRAM 
Cc 

P4=PI/180 

P5=P1¥*PG 

P6=PZ*PG 

P7=P3¥PG 

S=1 
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DO 200 I=0,90 

TH=I*P4¢ 

DO 100 PHI=0,180 
PH=PHI¥PG 
W1=K*¥D1*SIN( PH J¥SIN( TH ) 

— — W2SK¥D2*SIN( PH J¥SIN( TH J-P5 

~ W3=K*D3*SIN( TH )*COS( PH )-P6 
WG=K¥DG¥SIN( TH )*¥COS( PH )+P7 
W=COS(W1 ) ¢A¥COS(W2 )+B*COS(W3 )+C¥COS( WG ) 
G=SIN(W1 J-AXSIN(W2 J)-B¥SIN(W3 )4C¥SINCWG ) 
H=SQRT(WXN+GXG )/(1+A+B+C ) 
2( PHI )=20*ALOG10(H ) 
IF((Z(PHI) .GT.-10. ).AND.(Z(PHI).LT.-9.)) THEN 

Z(PHIJ= -30. 


ELSE 
Z(PHIJ= 2Z( PHI) 
END IF 
2M=Z( PHI ) 
Cc WRITE(8,1) I,PHI,2M 


IF (2M .EQ. -30.) THEN 
V(S) = FLOAT (PHT) 
U(S) = FLOAT (T) 


C WRITE(8,1) V(S),U(S),2M 
S=S+l1 

END IF 
Cc WRITE(8,1) I,PHI>Z 
Cl FORMAT( 2X%,2(13,5X),F10.5) 
C WRITE(8,1) AX(L),BY(L),2M 
Cl FORMAT (2X%53(F10.355X)) 
100 CONTINUE 
200 CONTINUE 
Cc 
C www ee ew emer eww eee ee ee ee ee een come eee wee’ coe e eee eee er ee eo cee cceeeceeececcca 
A PROGRAM FOR PLOTTING BY SHERPA OR BY TEK618 
C 


CALL MEDBUF 
C CALL COMPRS 
CALL TEK618 
CALL PHYSOR(1.,1.) 
CALL PAGE(15.5,12.) 
CALL AREA2D(10.5,9.1) 
CALL FRAME 
CALL COMPLX 
CALL HEIGHT( .2) 
CALL XNAME('PHI(DEG) $',10) 
CALL YNAME( ‘THETA(DEG)$',10) 
CALL YTICKS(2) 
CALL XTICKS(4) 
CALL GRAF(0,20,180,0,10,90) 
CALL CURVE(YV Rigesat) 
CALL MESSAGI 'THETA=25 DEGREE',15,11.,8.5) 
CALL MESSAG( 'K=5.5',5,11. 58.2) 
CALL MESSAG( 'D1l=1. D2=1.',11,11.57.9) 
CALL MESSAG( 'D3=1. DG=1.',11,11.57.6) 
CALL MESSAG( 'P2=180 P3=0 PG=180',18,11.,7.3) 
CALL MESSAG( 'A=1. B=l. C=1.'5,14511.57.0) 
CALL RESET( "HEIGHT ) 
CALL HEIGHT(.4) 
CALL MESSAG( 'TRACK( -20.0B?',13,3.1,9.8) 
CALL RESET('ALL' ) 
CALL NOCHEK 
CALL ENDGR(0O ) 
CALL PHYSOR(11.85,7.6) 
CALL AREA2D(3.5,2.5) 
CALL FRAME 
CALL ENDPL(0) 
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APPENDIX E 





P3D 
* * 
% THE OBJECTIVE OF THIS PROGRAM IS TO CALCULATE AND DRAW % 
% THE DIRECTIVITY PATTERN IN THREE-DIMENSION. * 
* haa abe Gla METHOD IS THE SUBROUTINE SURMAT OF DISSPLA x 
* PLANE. * 
x x 





Z : DIRECTIVITY IN dB 

2MM : DOUBLE MATRIX TO GET 9-@ plane 

X : DUMMY VARIABLE OF THETA 

Y : DUMMY VARIABLE OF PHI 

2M : DUMMY VARIABLE OF DIRECTIVITY 

Pl : PHASE (DEGREE) OF SOURCE A 

P2 : PHASE (DEGREE) OF SOURCE B 

PZ : PHASE (DEGREE) OF SOURCE C 

PG : VARIABLE TO CHANGE DEGREE TO RADIAN 

PS : RADIAL VALUE OF Pl 

P6 : RADIAL VALUE OF P2 

P7 : RADIAL VALUE OF P2 

K : WAVE NUMBER 

Dl : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE 2 
D2 : DISTANCE: FROM THE GEOMETRICAL CENTER TO THE SOURCE A 
D3 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE 8B 
D4 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE C 
PHI : PHI IN DEGREE 

I : THETA IN DEGREE 

PH : PHI IN RADIAN 

TH : THETA IN RADIAN 

PI : 3.1616917 


JBIBRBEBHRRRBBRBEEHRHHHE VARIABLE DECLARATION 288GQRRBRBBEHBEBRERHHEHEHRHE 
DIMENSION Z( 2500) ,ZMM( 46,46) 
REAL X»sY»ZM»P4,P5,P6,P7,W1 »W2»W3 »>W4,WN,G)H»K,01,02,03,04,PH,TH,PI 
real a»byc 
INTEGER PHI,P1,P2,P3,I 


KKK KKK KK KK KK KK KK KOK OK OK OK OK OK OX 
KKK KKK KKK KKK KK KKK KKK Ky KK 


c 
C em wm mew ew eww enw wow ew wo wo ee ww wo ee ww ww ww wn ww ow ee eee ewe 
: A PROGRAM FOR PLOTTING BY SHERPAOR BY TEK618 
Cc 
CALL COMPRS 
CALL TEK618: MOVE C TO THE UPPER LINE WHEN DRAWING WITH TEK618 
CALL RESET (3HALL) 
CALL PHYSOR(0.5>1.) 
CALL PAGE(15.5,12. ) 
CALL SCMPLX 
CALL AREA2D(10.5,9.15) 
CALL X3NAME( "THETA( DEGREE )',13) 
CALL Y3NAME( ‘PHI( DEGREE )’,11) 
CALL Z3NAME( 'DIRECTIVITY(0B)',15) 
CALL VOUEMSD (17.517. 517.) 
CALL VUANGL( +20. ,+20.,500. ) 
CALL INTAXS 
CALL ZAXANG (-90. ) 
CALL GRAF30(0.,5.,90.,180. »6.,360. »-40. »2.,0. ) 
CALL BLSUR 
? CALL BGNMAT( 46,46 ) 
C ww mw eee eww mmm m wom ome ewww www www wm ewww’ wm em em eee woe ee eee wee ee wees eee eee eee 
: INPUT PARAMETERS 
C 
PI=3 .1415917 
K=5.5 
Di=l. 
02=1. 
03=1. 
04=1. 
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P3= 180 
A = Hl e 
B=l1. 
c= il. 
c 
Cc se eee eee Seee eee ee ee eeee ee eee eee ee ee ee eee ee Geese ee eeeeeeeee eee eee es eesceeee 
c MAIN PROGRAM 
c ee Se See Seee Seee See eeeeeeee coecee es eee eee eee ee ee ee ee eee eee ee eee eeeoeacuse 
Cc 
PG=PI/180 
PS=P1¥P4 
P6=P2*PG 
P7=P3¥PG 
DO 200 I=0,90;2 
TH=I¥PG 
DO 100 PHI=180,360 »4 
PH=PHI*¥PG 
W1L=K¥D1L*¥SIN( PH J¥XSIN( TH ) 
W2=K*¥D0 2*SIN( PH )¥SIN( TH )-PS5 
W3=K*¥D3*SIN( TH )¥COS( PH )-P6 
WG=K¥DGXSIN( TH )¥COS( PH )4+P7 
W=COS(W1 )+A¥COS(W2 )+4B¥COS(WS3 )+CXCOS( WG) 
G=SIN(W1 J-AXSIN(W2 J-B¥XSIN(W3 J 4CXSINOWG ) 
H=SQRT (W¥WN4+GXG )7(14A4B4C ) 
Z( PHI )=20*ALOG1O(H ) 
IF(Z( PHI ).LE.-¢0.) THEN 
Z( PHI )=-40. 
LSE 
Z( PHI J= Z(PHTI ) 
ENO IF 
x=I 
Y=PHI 
Z2M=Z( PHI ) 
WRITE( 8,1) X»Y»52M 
1 FORMAT(2X,3(E10.3,5X)) 
CALL GETMAT(X,Y,2Z2M,1,0} 
100 CONTINUE 
200 CONTINUE 
c 
c S82 OS Se G22 eS f OSS SSBF See se SF 2S 22 SS e222 SBF S28 2 TSF ee2e2 Sec eF ee e2e eee eee 2 eee eee 
: A PROGRAM FOR PLOTTING BY SHERPAOR BY TEK618 
c 


CALL ENOMAT( 2ZMM,0) 
CALL SURMAT (2MM,1,465,1,46,0) 
CALL DASH 
CALL MARKER( 17) 
CALL COMPLX 
CALL HEIGHT(.2) 
CALL SURVIS( ‘BOTTOM’ ) 
CALL MESSAG( 'K=5.5',5,11.,8.5) 
CALL MESSAG( 'O1=1. O2=1.',11,11. ,8.2) 
CALL MESSAG( 'D3Z=1. OG=1.',11,11.,7.9) 
CALL MESSAG( ‘VIEW ANGLE= 20, 20',18,11.,7. 
CALL MESSAG( 'P2=180 P3=0 PG=180',18,11.,7. 
CALL MESSAG( °A=1. B=1. C=1.',1¢,11.,7.0) 
CALL RESET( ‘HEIGHT' ) 
CALL HEIGHT(.4) 
CALL MESSAG( 'DIRECTIVITY PATTERN' ,19,3.05,9.8) 
CALL ENDGR(0O) 
CALL PHYSOR(11.35,7.6) 
CALL AREA2D0(3.5,2.5) 
CALL FRAME 
CALL DONEPL 
CHECKSUM=012362 
STOP 
ENO 


6) 
3) 
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APPENDIX F 





TRACE 

JEEBEREREREEERRHHHEHHEeHEE «=INFORMATION 

*% * 
* THE OBJECTIVE OF THIS PROGRAM IS TO CALCULATE AND DRAW % 
* THE CONSTANT dB LEVEL OF THE THREE-DIMENSIONAL PLOT. x 
* THE DRAWING METHOD IS THE SUBROUTINE CURVE OF DISSPLA. x 
¥ *% 
IBHEEHEEREESEHHEHReRReE VARIABLE DEFINITION 28888R8RRBHHRBORBREHREHHHHHHE 
x x 
* 2 : DIRECTIVITY IN p * 
* v : ARRAY TO TRACK x 
* u;: ARRAY TO TRACK * 
* Pl : PHASE (DEGREE) OF SOURCE A x 
% P2 : PHASE (DEGREE) OF SOURCE B * 
* P3 : PHASE (DEGREE) OF SOURCE C * 
% P4 : VARIABLE TO CHANGE DEGREE TO RADIAN * 
* P5 : RADIAL VALUE OF Pl cs 
% P6 : RADIAL VALUE OF P2 * 
* P7 : RADIAL VALUE OF P3 * 
* K : WAVE NUMBER x 
* D1 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE 1 % 
* D2 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE A * 
* D3 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE B % 
* DG : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE C % 
* PHI : PHI IN DEGREE % 
* I: THETA IN DEGREE * 
* PH : PHI IN RADIAN * 
* TH : THETA IN RADIAN * 
% PI s 3.1615917 x 
*% * 


DIMENSION 2(0:2500 ),ZMM(46,46),V{2000 ) ,U{ 2000 3 

REAL Po5T »sRoXaV¥ »>ZMyZZMoA5BoC oV5U 

REAL P4,P5,P6,P7>,W1 WZ WN3 »WNG»WN,G,H,K,D1,02,D03 D4,PH,TH,PI 
INTEGER PHI,P1,P2,P3,I>,S 


Cc 
Cone ee a a aes ee a ee oe ee eee a ee ce 
C INPUT PARAMETERS 
6 aaa eae aa ae a a ae a ee ee ee ee 
C 
PT=3.1415917 
K=10.55 
D1=0.5 
02=0.5 
03=0.5 
D4=0.5 
Pl= 0 
P2=180 
P3=180 
= Ls 
B=l. 
c=l. 
C 
(Coon n sees we mescaswcnccescee see eee oe ee- asa e eee eae eee ee 
. MAIN PROGRAM 
C 
P4=PTI/180 
P5=P1*PG 
P6=P2*PS 
P7=P3*P4 
S=1 
DO 200 I=0,90 
TH=I¥PS 
DO 100 PHI=0,180 
PH=PHI*P4 


W1=K¥DI*SIN( PH J¥SIN(TH ) 
W2=K¥DZ*¥SIN( PH JXSIN( TH J-PS 
W3=K*¥DS*SIN( TH )*COS( PH )-P6 

WG =K*¥DG*SIN( TH )*COS( PH 3+P7 

W=COS(W1 )+A®COS(W2 )+B¥COS(W3 )+C¥COS( WG ) 
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GSSIN(W1 )-AXSIN(W2 )-BXSIN(W3 )+C¥SINI WG ) 

H=SQRT ( WXN+GXG )/7(14+A+B4C ) 

Z( PHI )=20%ALOG10(H } 

IF((Z(PHI) .GT.-10.).AND.(Z2(PHI).LT.-9.)) THEN 


Z(PHIJ= -30. 
ELSE 
Z(PHI)= Z(PHI) 
END IF 
2M=2( PHI ) 
C WRITE(8,1) I,PHI,2M 


IF (2M .EQ. -30.) THEN 
V(S) = FLOAT (PHI) 
U(S) = FLOAT (TI) 


Cc WRITE(8,1) V(S),U(S),2M 
S=S4+1 
END IF 
C WRITE(8,1) I,PHI,2 
Cl FORMAT (2X ,2(13,5X),F10.5) 
Cc WRITE(8,1) AX(L),BY(L),2M 
Cl FORMAT ( 2X»3(F10.3,5X)) 
100 CONTINUE 
200 CONTINUE 
f 
. A PROGRAM FOR PLOTTING BY SHERPA OR BY TEK618 
Cc 


CALL MEDBUF 

CALL COMPRS;MOVE C TO TEK619 TO PRINT USING SHERPA 
CALL TEK618 

CALL PHYSOR(1.,1.) 

CALL PAGE(15.5,12.) 

CALL AREA2DE10.559.1) 

CALL FRAME 

CALL COMPLX 

CALL HEIGHT(.2) 

CALL XNAME('PHI(DEG) $',10) 

CALL YNAME( 'THETA( DEG )$',10) 

CALL YTICKS(2) 

CALL XTICKS(G) 

CALL GRAF(0,20,180,0,10,90) 

CALL CURVE(V,U;S5-1) 

CALL MESSAG( 'THETA=25 DEGREE',15,11.,8.5) 
CALL MESSAG( 'K=5.5',5,11. 58.2) 

CALL MESSAG( 'D1=1. D2=1.',11,11.,7.9) 
CALL MESSAG( 'D3S=1. D4=1.',11,11.>7.6) 
CALL MESSAG( 'P2=180 P3=0 PG=180',18,11.,7.3) 
CALL MESSAG( 'A=1. Bel. C=1.',145,11. 57.0) 
CALL RESET( ‘HEIGHT ' ) 

CALL HEIGHT( .4) 

CALL MESSAG( 'TRACK(-20.0B)',13,3.1,9.8) 
CALL RESET( ‘ALL' ) 

CALL NOCHEK 

CALL ENDGR(O) 

CALL PHYSOR(11.85,7.6) 

CALL AREA2D(3.5,2.5) 

CALL FRAME 

CALL ENDPL(0) 
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APPENDIX G 
LIST OF COMPUTER OUTPUT 


@ SAMPLE CASE 1 


F999 IEEE EE EEE III EIEIO IEE 
A B eC K D D, oO D. D, @ ©, oO 6 “ F 


1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 Oo 180 180 0 31 0.955567 -0.394776 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 52 0.970322 -0.261679 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 33 0.982069 -0.157164¢ 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 34 0.990826 -0.080049 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 35 0.996641 -0.029226 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 36 0.999579 -0.003656 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 O 180 180 0 37 0.999728 -0.002363 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 38 0.997193 -0.024419 
1.0 1.0 1.0 10.55 0.50 0.50 9.50 0.50 0 180 180 0 39 0.992094 -0.068947 
1.90 1.0 1.0 10.55 0.50 0.50 0.50 0.50 O 180 180 0 40 0.984566 -0.135104 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 41 0.974756 -0.222079 
1.0 21.0 1.0 10.55 0.50 0.50 0.50 0.50 O 180 180 0 42 0.962821 -0.329093 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 O 180 180 0 41 0.948922 -0.455388 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 oS 0.933230 -0.600225 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 45 0.915917 ~0.762881 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 O 180 180 0 46 0.897156 -0.942645 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 O 180 180 0 47 0.877121 -1.138807 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 48 0.855985 -1.350672 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 49 0.833918 -1.5775355 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 50 0.811083 -1.818697 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 81 0.261312 -11.656804¢ 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 O 180 180 0 82 0.255355 -11.857100 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 O 180 180 0 83 0.250133 -12.036570 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 84 0.245633 -12.194245 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 85 0.241845 -12.329250 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 86 0.238759 -12.403822 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 87 0.236367 -12.528270 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 88 0.234663 -12.591125 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 89 0.233642 ~12.628973 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 6180 6 6180 0 90 0.233302 -12.641622 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 6 180 180 0 aL 0.233642 -12.628973 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 92 0.234663 -12.591125 
1.0 1.0 1.0 10.55 0.50 0.50 90.50 0.50 O 180 180 0 93 0.236367 -12.528270 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 94 0.238759 -12.440822 - 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 95 0.241845 -12.329267 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 96 0.245633 -12.194263 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 97 0.250133 -12.036578 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0.6180 «6180 0 98 0.255355 -11.857116 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 he) 0.261311 -11.656837 
1.0 1.0 1.0 10.55 0.50 0.50 0.50 0.50 0 180 180 0 100 0.268014 -11.436835 
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